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Resume 

The microscopic nature of an insulating state in the vicinity of a superconducting state, in the presence 
of disorder, is a hotly debated question. While the simplest scenario proposes that Coulomb interactions 
destroy the Cooper pairs at the transition, leading to localization of single electrons, an alternate possibility 
supported by experimental observations suggests that Cooper pairs instead directly localize. The question 
of the homogeneity, granularity, or possibly glassiness of the material on the verge of this transition is 
intimately related to this fundamental issue. Here, by combining macroscopic and nano-scale studies of su¬ 
perconducting ultrathin NbN films, we reveal nanoscopic electronic inhomogeneities that emerge when the 
film fhickness is reduced. In addifion, while fhicker films display a purely fwo-dimensional behaviour in fhe 
superconducfing flucfuafions, we demonsfrafe a zero-dimensional regime for fhe fhinner samples precisely 
on fhe scale of fhe inhomogeneifies. Such behavior is somehow infermediafe befween fhe Fermi and Bose 
insulafor paradigms and calls for furfher invesfigafion fo undersfand fhe way Cooper pairs confinuously 
evolve from a bound sfafe of fermionic objecfs info localized bosonic enfifies. 
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Understanding the mieroseopie proeesses oeeurring at the supereonduetor-insulator transition 
(SIT) in ultrathin films is important not only for fundamental reasons, but also for applieative 
purposes mm. The mierostruetural properties are known to play a key role and the samples ean 
then be divided into two groups Q : granular thin films and homogeneously disordered thin films. 
For the former, the SIT is driven by the eompetition between the intergrain Josephson eoupling, 
favoring the deloealization of pairs, and the Coulomb eharging energy of the grains, whieh renders 
eharge fluetuations energetieally expensive IH. In the ease of nominally homogeneously disor¬ 
dered films however, whieh are the objeet of this Artiele, several seenarios have been proposed. 
On the one hand, what is often referred to as the "fermionie" seenario proposes that the meeha- 
nism that drives the transition is the redueed sereening of the Coulomb repulsion with inereasing 
disorder, weakening pairing and redueing the eritieal temperature Tc [01, as observed in [j6l. In 
this ease, the insulating state is eomposed of loealized fermions and, in partieular, standard para- 
eonduetive fluetuations are expeeted due to Gaussian-distributed short-lived Cooper pairs. On the 
other hand, pairs may survive the SIT in a "bosonie" seenario, in whieh the gap persists above Tq 
despite the loss of phase eoherenee. In this framework, the bosonie pairs either loealize beeause 
disorder-enhaneed Coulomb interaetions destroy their phase-eoherent motion at large seales [|7l[8]| 
or disorder itself ean blur the pair phase eoherenee without any relevant role of the Coulomb re¬ 
pulsion [l9l-[T^ . On the experimental side, it has been shown through eareful study of Little-Parks 
oseillations, that either fermionie [[T3ll or bosonie [|T4lllT5l transitions may oeeur. In the latter ease, 
it was also proposed that the supereondueting state is eharaeterized by an emergent disordered 
glassy phase [[T^ with filamentary supereondueting eurrents [|9||. An anomalous distribution of the 
supereondueting order parameter was proposed by theorists [TT^fT^ . and observed experimentally 
ifTTl ITSi . A numerieal approaeh to uniformly disordered supereonduetors IIT^ has also suggested 
that there is a eontinuous evolution |j20l from the weak disorder limit, were the system has a rather 
homogeneous fermionie eharaeter, to the strong disorder limit, where marked inhomogeneities ap¬ 
pear in the supereondueting order parameter, with an emergent bosonie nature eharaeterized by 
a single-partiele gap persisting on the insulating side. A great deal of experimental aetivity has 
been devoted to the more disordered part of the SIT IfTTl 1211 while the intermediate region where 
fermionie Cooper pairs begin to evolve into bosonie pairs has not been extensively aeeessed. The 
aim of this work is preeisely to fill this gap by reporting experiments whieh shed light on how 
Cooper pairs evolve with inereasing disorder, giving rise to ineipient inhomogeneous bosonie fea¬ 
tures. More speeifieally, we present here a study on a set of ultrathin NbN films that are nominally 
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homogeneous, but where eleetronie inhomogeneities and pseudogap emerge as the thiekness is 
redueed, together with experimental indieations in favor of fermionie meehanisms. Indeed, while 
the thiekest films {d > 2.2 nm, Tq > 0.3 are found to be rather homogeneous ^22} with 
two-dimensional (2D) Aslamazov-Larkin (AL) supereondueting fluetuations, the thinnest samples 
offer a more eomplex behavior eharaeterized by inhomogeneous supereonduetivity, indieation for 
a pseudogap above Tc (a seemingly "bosonie" feature), in agreement with the literature ifTTlfTSlI^ 
and by the establishment of a zero-dimensional (OD) regime of gaussian supereondueting fluetua¬ 
tions (a fermion-like hallmark). We show that for films below this threshold film thiekness and 
eritieal temperature, the supereondueting fluetuations measured by transport above Tc behave in 
agreement with a formal OD limit of the AL theory of paraeonduetivity, in a substantial range of 
redueed temperature e = log(T/Tc). This indieates that these fluetuations are still eonventional 
and seemingly BCS-like, but somehow eonfined in a "supergrain" over a lengthseale 1. The super¬ 
eondueting inhomogeneities as evideneed by seanning tunneling speetroseopy at low temperature 
eorrespond to eleetronie domains of typieal size Lj/2 preeisely of the order of I ~ 50 nm i.e., 
mueh larger than any definite struetural seale of the system. The paradoxieal presenee of AL flue¬ 
tuations together with indieations for bosonie features sueh as the pseudogap, well established in 
these systems ifTTlfT^I^ . is one of the most intriguing results of this work. We suggest that these 
two features ean indeed be reeoneiled if the pseudogap in our system has a fluetuational origin 
[|24ll . We propose that the pronouneed amplitude of the pseudogap observed in the ultrathin films 
that are the objeet of our investigation arises from a slowing down in the diffusion of the fluetua- 
ting Cooper pairs, whieh exhibit a tendeney to loealize on the typieal inhomogeneity seale Lj and, 
at the same time, give rise to the OD AL behavior. 

Results 

Our samples eonsist of ultrathin NbN films grown ex situ on sapphire substrates. Details of the 
fabrieation proeess may be found in [|25ll . (See also the Appendix). The different samples studied 
together with their thiekness, eritieal temperature, and resistanee per square at room temperature 
may be found in Table 1 of the Appendix. 

Transmission Eleetron Mieroseopy (TEM) measurements of our NbN films (see Eig.5 in the 
Appendix) were able to show that the films are smooth and well erystallized, and ean be viewed 
as a elosely paeked assembly of nanoerystallites of different orientations. These eontiguous na- 
noerystals have typieal lateral dimensions dg of about 2—5 nm. The eharaeteristie film strueture 
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Figure 1. Local inhomogeneities of the superconducting state, (a) Topographic map of the NbN area under 
study for sample Xq (Tq = 3.8 K). (b) Corresponding map displaying the superconducting gap inhomo¬ 
geneities at 300 mK. Gap inhomogeneities appear on a much larger scale than the size of the nanocrystals 
{dg ~ 2 — 5nm) constituting the NbN films and seen in the topographic map.(c ) Cross-correlation map 
between the topographic and spectroscopic maps revealing the absence of correlation, (d) Autocorrelation 
map of the gap map (b). (e) Schematics of the relevant length scales, (f) Radial profile exfracfed from fhe 
aufocorrelafion map (d). The correlafion lengfh L* is defined as fhe abscissa of fhe firsf peak away from fhe 
cenfer and so is abouf lOOnm in fhe presenf sample. The cenfral peak widfh, abouf Lj/2 ~ 50 nm gives an 
esfimafe of fhe fypical domain size. 

is also reflected in topographic scanning tunneling microscopy images, as the one presented in 
Fig.[^ for a nominally 2.14nm-thick sample Xq (Tc= 3.8 K). The film surface is very flat and the 
observed nanoscale structures correlate well with the nanocrystals revealed by TEM. At the same 
time the landscape of the sample also displays smooth inhomogeneities on a larger scale of a few 
tens of nanometers. 

Probing the inhomogeneities with STM and STS — In order to get insight into the superconduc¬ 
ting (in)homogeneity of these ultrathin films, we performed scanning tunneling microscopy (STM) 
and spectroscopy (STS) experiments. Typical results are displayed in Fig.[2 Fig.[2a) shows a to- 
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pographic map of the area under study for sample Xq (Tc = 3.8 K). In Fig. [^b) we report the 
extraeted supereondueting gap map A(X, Y) at 300 mK i.e. well inside the supereondueting state 
for the area eorresponding to the topographie map of Fig. [^a) . A (A, y) is defined from the dl/dV 
speetra as the peak-to-peak energy (see Fig. 6 in the Appendix for the eorresponding representa¬ 
tive dl/dV speetra). In agreement with previous measurements in similar systems ifTTl 1^ . we 
observe spatial variations of the gap (see Fig.[2b)) and of the integrated in-gap eonduetanee (see 
Appendix). However, we reveal here for the first time that these gap variations are not correlated 
to the topography of the surface. Indeed, the eross-eorrelation map between (a) and (b) reveals the 
absenee of eorrelation [see Fig.[2e)]. 

In Fig. [^d) we report the autoeorrelation map of the gap map shown in Fig. [^b). It allows 
us to extraet both the typieal size of domains of eonstant gap values, hereafter ealled supergrains, 
and the typieal distanee Lj between the eenters of sueh adjaeent supergrains. The radial profile 
extraeted from the eenter of the autoeorrelation map is shown in Fig.[2f), featuring a eorrelation 
length Li of about 100 nm and a typieal domain size of about Tj/2. This is more than ten times 
larger than the size of the mierostruetural grains dg ~ 2 — 5 nm (see Fig.l^e) depieting the relevant 
lengths seales). Above the global eritieal temperature Tc these inhomogeneities persist, as well as a 
gap-like feature ealled pseudogap. In addition, our measurements reveal that the energy integrated 
in-gap eonduetanee maps at 300 mK and at 4.2 K, performed on the same area of Fig. [^a) are 
eorrelated (see Fig. 6 in the Appendix). We deduee that at 4.2 K Lj remains roughly the same, 
whilst the speetroseopie eontrast between the regions ehanges with temperature. 

By eontrast, similar speetroseopie studies earned out on thieker samples (d > 2.3 nm, Tc > 
0.3 revealed a mueh more homogeneous supereondueting phase [|22l . Therefore, some short- 
seale inhomogeneity in the supereonduetive properties emerges for the thinner samples at a seale 
Li while the system is strueturally homogeneous over the same length seale. This is eonsistent 
with predietions from Monte-Carlo simulations |[T9l . 

Analysis of the paraconductivity — In order to probe the influenee of these nanoseale inhomo¬ 
geneities on the supereondueting thermal fluetuations, we performed transport measurements in 
the vieinity of Tc and extraeted the paraeonduetanee per square Xa{T) = a{T) — crAr(T), i.e., 
the exeess eonduetanee per square due to supereondueting fluetuations in the normal state. Here, 
cr(T) is the square eonduetanee measured under zero magnetie field and (TAr(T) is the normal state 
square eonduetanee. The resistanee per square is displayed on Figj^for the different samples, toge¬ 
ther with the extrapolated [(a) solid lines] or measured [(b) full symbols] resistanee of the normal 
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Figure 2. (a) Square resistance under 0 T as a function of temperature for the thicker samples B2 {Tc = 
7.1 K, +), Cl {Tc = 9.4 K, x), Fq {Tc = 9.0 K, *), with the corresponding extrapolated normal state 
resistance (solid lines, see text), (b) Square resistance under 0 T as a function of temperature for the thinner 
samples Yq {Tc = 4.3 K, open circles), Xq {Tc = 3.8 K, open squares), A 2 {Tc = 4.5 K, open diamonds), 
and Ai (aged A 2 , Tc = 2.4 K, open triangles). The square resistances under a perpendicular magnetic field 
of 14 T are reported with the corresponding full symbols. The solid lines are for the extrapolated normal 
state resistances. 

state. (See the methods Seetion for details on the determination of the normal state.) 

In Fig.j^a) the variation of Act with the reduced temperature e = \ n { T / Tc ) is shown for 
samples B 2 , Ci, and Fq : the observed critical exponent —1 is consistent with the Aslamasov- 
Larkin prediction for 2D systems (AL 2D) [l26ll Act = e^/ { 16 he ) in the range 0.02 ^ e ^ 0.2 —0.9, 
as it was reported previously for these films ||25l. Remarkably, the extracted experimental AL 
2D prefactor matches the theoretical one, without any adjustable parameter. This suggests that 
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the fluctuations in this case are BCS-like, and that Maki-Thomson (MT) fluctuations [l27l 1^ or 
density-of-state (DOS) corrections [ffl |29l [301 are absent or negligible. 

Proceeding in a similar way, we extracted the paraconductivity of the thinner samples A 2 , A 4 
(aged A 2 ), Xo and Vq (shown as open symbols in Fig.j^b). The AL 2D prediction is also displayed 
in the same Figure (purple solid line). The experimental paraconductance is found to deviate stron¬ 
gly from the AL 2D behavior over a significant temperature range, even when using Tc as a free 
parameter, and to exhibit over a substantial range of reduced temperature a very specific law. 
Act ~ corresponding to OD fluctuations, previously observed in granular materials [|3TI - l34l . 
An empirical fitting function. Act = 0.03e^/ (he^) for A 2 and Xq is plotted for comparison (ma¬ 
roon solid line). Similar behaviors, i.e., Acr = 0.02e^/(^e^) for Yq, and Act = 0.015e^/(^e^) for 
A 4 , were found to hold for the other samples. The paraconductance data, e.g., for sample A 2 also 
suggests a 0D-2D crossover [see also the lower inset in Fig.|^b)] as previously observed in Il35ll 
and discussed in [|^ . We point out that the fluctuational critical temperature may differ in the 
two regimes, so that the crossover is more evident when the pararesistivity Xa~^ is plotted as a 
function of T, without making any choice for Tc, rather than e, which depends on Tc [see the 
upper inset in Fig.j^b)]. 

Concerning the absence of MT terms, we stress out that, with pair breaking arising only from 
electron-electron interactions, MT paraconductance is less singular than the AL term in 2D [|^ . 
and even less so for OD AL. On the other hand, the presence of a sizable pseudogap suggests that 
DOS corrections should be present. DOS corrections, however, lead to a decrease of the paracon¬ 
ductance. In our case, instead, the paraconductance in the pseudogap regime is found to be even 
more singular with dependence that can not be explained by DOS contribution. This clearly 
indicates that DOS corrections, although expected, are immaterial in this case. 

Study of the magnetic field driven transition — The above study near Tc was complemented by 
the analysis of the transport properties of the thinnest samples at the transition to the normal state 
driven by magnetic field, well below Tc- We found that the curves R{H) cross at a specific point 
[Re] He} (see also Fig. 8 in the Appendix for further details). Following Ref. [371 we analyzed 
the curves in the vicinity of this point and evidenced a scaling behavior of the type R/Rc = 
f{\H — Hc\T~°') with the critical exponent a ^?>/2 (See the data in Fig.|^a) for sample A 2 ). The 
occurrence of such scaling behavior, where T is the only relevant scale, marks the existence of a 
quantum critical point (QCP) at zero temperature and for H = He- Consequently, the exponent 
a can be expressed as a = l/{uz), where u is the exponent that rules the variation of the spatial 
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Reduced Temperature 


Reduced Temperature 


Figure 3. (a) Extracted square paraconductance for the thicker samples B 2 (Tc = 7.1 K, +), Ci (Tq = 
9.4 K, x), and Fq {Tc = 9.0 K, *) as a function of the reduced temperature e = ln(r/Tc). The agreement 
with the Aslamasov-Larkin prediction for a 2D system (purple solid line) is excellent, without any adjustable 
parameter, (b) Extracted square paraconductance for the thinner samples Yq {Tc = 4.3 K, open circles), Xq 
{T c = 3.Sir, open squares), A 2 {Tc = 4.5 K, open diamonds), A 4 {Tc = 2AK, open triangles). The 
maroon solid line corresponds to Act = 0.03e^/ (he^). The expected AE 2D square paraconductance is also 
shown (purple solid line). The upper inset shows the pararesistance as a function of temperature for 

the sample A 2 , emphasizing the 2D [oc (T — Tc)] and OD [oc (T — Tc)^] behaviours, the crossover form 
one to the other, as well as the different fluctuative critical temperatures in the two regimes (T^^ = 4.9 K, 
T^^ = 4.4 K). The lower inset shows the crossover to the 2D AE square paraconductance for the sample 
A 2 as a function of the reduced temperature. (In this case e = ln{TjT^P) was used.) 
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correlation length ^ ~ |if — and z is the dynamieal eritieal exponent ~ 1/T. 

In order to extraet the values of v and 2 ; from a ~ 3/2, similar measurements were performed 
at a fixed temperature, T =1.9 K, for different values of the eleetrie field, i.e., of the bias voltage 
aeross the sample [see the data for sample in Fig. |^b)]. Here again, the R[B.) eurves exhibit 
a eommon erossing point eorresponding to the QCP (see Fig. 9 in the Appendix). The sealing 
analysis in the vieinity of this point, with a sealing funetion of the form R/Rc = g{\H —Hc\E~^), 
yielded (3 ^ 3/4. Following the analysis in Refs. lT7l - [39l we express (3 as (3 = l/i'iz + 1). The 
two independent determinations of uz and u{z + 1) allow us to establish u = 2/3 and z = 1. The 
latter is preeisely the result expeeted, e.g., in systems with (weakly sereened) long-range Coulomb 
interaetions, while the former is eonsistent with a 3D elassieal XY model or 2D quantum XY 
model (with z=l). A similar analysis, performed on sample Fq under high pulsed magnetie field 
allowed to evidenee a plateau between 1.5K and 8K for He = 18.6K and to extraet a produet of 
eritieal exponents uz ~ 2/3. 

Discussion 

As we report in the Appendix, the eoherenee length exponent that leads to the observation of 
an anomalous power law Aa ~ e“^, is eonsistent with a formal ealeulation of AL fluetuations 
in a OD system. Act is eonverted into the measured paraeonduetanee per square by means of 
a suitable length seale /, whieh represents the size of the OD fluetuating domains, in the plane 
parallel the film, yielding 5aD=o = Dedueing for ,^0 a value of 5.5 ±0.5nm from the 

value of He at the QCP, in agreement with extrapolated values in [[251, it is possible to extraet I 
from the paraeonduetivity data. One finds, e.g., I = 28 nm for samples A 2 (Te = 4.5 K) and Xq 
(Te = 3.8 K), I = 35 nm for sample Yq (Te = 4.3 K), and I = 40 nm for sample A4. This length I 
is in good quantitative agreement with the typieal domain size Lj/2 ~ 50 nm extraeted from STS 
data at 300mK and at 4.2 K for sample Xq (Te = 3.8 K). This means that the length Li/2 ~ I 
and not the real grain size dg <C I sets the seale for the OD fluetuating domains, until ^ beeomes 
larger than Lj. Further theoretieal investigation is neeessary to understand the physieal meaning of 
the presenee of these OD domains, that behave like smooth "supergrains". The observed behavior 
eould aetually be interpreted in terms of slower diffusion of the Cooper pairs thus inereasing their 
lifetime on length seales smaller than 1. By eontrast, a standard diffusive behavior is reeovered 
at longer times and larger distanees giving way to standard 2D behavior eventually ruling the 
transition (see the Appendix for more details). 
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Figure 4. (a) Normalized resistance R/Rc of sample A2 as a function of the scaling variable \H — Hc\t 
for different temperatures. He = 10.7T and Rc = llkfl at the crossing point, t = was adjusted 

in order to obtain the best collapse of the data. Inset : Log-log plot of the parameter t vs temperature, used 
to determine the value oi vz = 2/3. (b) Normalized resistance R/Rc of sample A2 as a function of the 
scaling variable \H — Hc\u for different electrical fields, measured at T =1.9 K. u = was 

adjusted in order to obtain the best collapse of the data. Inset : Log-log plot of the parameter u vs. voltage 
used to determine the value o^v{z + l) = 4/3. 
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To our knowledge, the only previous evidence in the literature of a OD fluctuations regime in 
transport is for nominally granular or filamentary systems OTl [331 [34||. The novelty lies here in 
the observation of such behavior in a compound where the inhomogeneity arises in a “mild” way : 
The films are far from granularity because the OD behavior does not occur on the small scale of 
the crystallites, but rather on the larger typical scale I, comparable to the correlation length Li 12 
inferred from STS. The emergent (as opposed to the structural) character of the OD regions is 
also suggested by the lack of any correlation [see Fig. [Jc)] between the inhomogeneous domains 
observed with STS [Fig. [^b)] and the large-scale structural disorder observed in the topography 
of the crystallite ensemble [Fig.[2a)]. 

For the superconducting transition to be probed in transport, the OD fluctuations regime has 
to finally evolve to higher-dimensional behavior. Close enough to Tc a crossover to 2D behavior 
must (and does) occur, when the superconducting coherence length becomes larger than Lj and 
enables to couple different OD domains, following a scenario analogous to the Lawrence-Doniach 
description for lamellar materials [|^|40]| . As a matter of fact, the 0D-2D crossover is visible for 
sample Xq (d = 2.14nm), as well as for the sample A 2 {d = 2.16 nm). Please see the insets of Fig. 

However, we propose here a different physical explanation for this crossover. We start with an 
expression of the AL paraconductivity in D dimensions (see details in the Appendix) : 

Moil) 


daoie) = 


vre 

AhD 


d7 


(e + 7)s 


( 1 ) 


with a suitable "density of states" (weighted with current vertices) Afnil)- Setting D = 2, the 
standard AL result is found for A /2 ( 7 ) = 7 / 7 ^, corresponding to diffusion of fluctuating Cooper 
pairs in two dimensions. However, if 


■^ 2 ( 7 ) = 


7/71, 

7/7r, 


for 7 < 7 , 
for 7 > 7 , 


( 2 ) 


corresponding to a slowing-down of the diffusion of fluctuating Cooper pairs above a threshold 7 , 
a OD behavior 


Sa ~ 


e^7 

16 he^ 


(3) 


is found for e ^ 7 . A comparison with the formal extrapolation of AL fluctuations to D = 0 
yields 7 = 47r(^o/0^- For e <C 7 , the standard 2D AL paraconductivity is recovered. 

Finally, at even lower temperature (below the 2D AL transition temperature), the system might 
be governed by percolation physics or, alternately, by Berezinski-Kosterlitz-Thouless behavior. 
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This behavior, if any, should occur on a very restricted range of temperature, because our mea¬ 
surements always display the paraconductivity of standard gaussian fluctuations. If the bosonic 
scenario was applicable to our films, the transition should be ruled by dephasing of the parrs and 
paraconductivity should mirror the characteristics of the vortex fluctuations in the BKT transition 
m. So, at first sight, our findings seem at odds with the existence of a pseudogap on the normal 
side of the transition as is well established in the literature and can be seen in the conductance maps 
and dl/dV spectra in the supplementary material, and which is usually ascribed to the strong loca¬ 
lization of Cooper pairs [|T^[23]I . Our work suggests a different explanation. The diffusion slowing 
down of the fluctuating Cooper pair in the "supergrains" increases their lifetime. This might occur 
either because they locally find a more suitable environment, like, e.g., a locally higher Tc, or, 
conversely, because they display an increased tendency to localize for those wavevectors corres¬ 
ponding to the largest inhomogeneities. This second possibility seems more likely because our fits 
suggest that the OD critical temperature (i.e. the “local” critical temperature) is slightly lower than 
the 2D large-scale Tq- The incipient localization and diffusion slowing down of the Cooper pairs 
seems to have a sizable effect on depressing the density of state at the Fermi level, thus suggesting 
the possibility of a substantial fluctuational pseudogap ll24ll [[29l . 

The magnetic-field driven transition at low temperature can be interpreted as magnetic field- 
induced dephasing of the OD supergrains thereby accounting for the critical exponents z/ = 2/3 
of an XY model in 2-1-1 dimensions, the additional "1" coming from the dynamical critical index 
z = 1. The value z = 1 is pertinent, e.g., to systems with (weakly screened) long-range Coulomb 
interactions [|42l. It is also consistent with numerical calculations based on a Boson-Hubbard mo¬ 
del [l43]|. Similar values for u were observed in Bi [[39l or NbSi [|4^ amorphous thin films, whereas 
a large number of studies of the SIT point towards a different universality class with z/ = 4/3 
Il3^ [451 - 1471 . a critical exponent consistent with classical percolation. Our findings are therefore 
in agreement with a picture of phase fluctuating OD supergrains coupled d la Josephson, for the 
magnetically driven transition. 

The question arises of the origin of the electronic inhomogeneities. Indeed, while the structural 
small scale inhomogeneity associated to the structural grains (2—5 nm) appears to be irrelevant, we 
deal with three larger scales over distances of tens of nanometers (Li ^ 100 nm): a) the electronic 
inhomogeneity of the pseudogap seen by STS [Fig. 1(b)]; b) the scale of the OD AL behavior 
seen in transport, and c) the topographic smooth landscape [Fig. 1(a)]. Although the scale b) is 
obtained from transport and cannot be easily connected to a spatial structure, it is quite tempting 
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to associate the eleetronie seales of pseudogap and OD transport to the seale on whieh Cooper 
pairs tend to loealize (before they eventually eondense on the infinite seale of the 2D transition). 
On the other hand, the faet that there is no eorrelation between the topographieal map and the 
supereondueting gap map proves that the strueture is not responsible for the gap inhomogeneities 
in a trivial way. i.e. it is not some loeal parameter variation (for instanee thiekness or stoiehiometry) 
that induees a loeally-eorrelated variation of the supereondueting properties. This does not mean 
that disorder or struetural inhomogeneities are irrelevant. They are most probably relevant, but in 
a eomplex manner sueh as, for instanee, disorder is relevant for loealization but loealization length 
is not simply the distanee between impurities. 

In any ease the main result of our work is that a OD physios is found to spontaneously emerge in 
transport measurements on the same length seales as the length seale for the supereondueting inho¬ 
mogeneities, muoh larger than the struetural orystallites, indioating that this struetural oomplexity 
has an eleetronie eounterpart produeing anomalous diffusion of the Cooper pairs. 

Conclusion 

The emergenoe of (glassy) inhomogeneous supereondueting phases out of homogeneously di¬ 
sordered films has been proposed both theoretieally [fT^ [T^ [T9]| and measured experimentally 
lfr^[2Tl[22ll . Our STS measurements have evideneed eleetronie inhomogeneities with typieal eor¬ 
relation length {Li ce 100 nm) and typieal domain size Lj/2 in the supereondueting state below Tc 
and in the supereondueting fluetuations above Tc in NbN ultrathin films. These inhomogeneities, 
whieh are mueh larger than the struetural grains of the film made of nanometer-sized nanoerys- 
tals, eorrespond to simultaneous spatial variations of the energy of the supereondueting gap and 
of the energy-integrated in-gap eonduetanee, and are found to beeome dominating for the thinner 
samples, where they strongly affeet the supereonduetive thermal fluetuations. Speeifieally, we have 
shown that, for films with nominal thiekness d < 2.3 nm or Tc < 0.3 these inhomogeneities 
are assoeiated with speeifie exponents in the dependenee of the gaussian supereondueting fluetua¬ 
tions with redueed temperature. These exponents are eonsistent with AL fluetuations eonfined into 
OD supergrains of size I ~ Tj/2. On the other hand, at low temperature, the analysis of the ma- 
gnetie field driven transition is eonsistent with a 2D quantum XY model, whieh is also evoeative 
of OD supergrains eoupled d la Josephson. In this regime, the Cooper pairs are tightly bound and 
the transition is ruled by supergrain dephasing, in eontrast to the transition at finite temperature 
where AL amplitude fluctuations are instead observed. 
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We propose that, above Tc, diffusion slowing down of the Cooper pairs in the supergrains 
inereases their lifetime and leads to a depression in the density of states at the Fermi level, yiel¬ 
ding a fluctuational pseudogap. The origin of these inhomogeneities as well as of this anomalous 
diffusion proeess is still to be investigated. In any ease, this work appeals for a new theoretical 
microscopic description for this peculiar state of matter, with interplay of localization and pairing, 
intermediate between the Fermi and the Bose insulators paradigms. 

Methods 

Scanning tunneling spectroscopy The electronic inhomogeneities in the superconducting state 
were probed by using scanning tunneling spectroscopy (STS) at 300 mK or 4.2 K and establishing 
full I{V) STS grids of typically 256x256 resolution over a 300x300 nm^ area. The dl/dV tunne¬ 
ling conductance spectra were obtained by numerical derivation of the I (V) data. The conductance 
maps shown in Fig. 6 of the Appendix were obtained by integrating in energy the individual dl/dV 
conductance spectra in an energy range of ±0.9 meV. The radial profile shown in Fig.[2f), extrac¬ 
ted from the autocorrelation map [Fig.[ 2 d)] of the gap map [Fig.[ 2 b)], is the circular average of all 
profiles passing through the center of the autocorrelation map. The measurements were performed 
with a Pt/Ir tip. 

Transport measurements Transport measurements were carried out by means of a standard four- 
probe technique in a 14 T Quantum Design PPMS (Physical Properties Measurement System). Re¬ 
sistance data was extracted using either a dc current or an ac current at 17 Hz. Four-probe contacts 
were made by directly bonding A1 wires to the NbN films. Six different samples were experimen¬ 
tally analyzed corresponding to different thicknesses : 2.10nm (Fq), 2.14nm (Xq), 2.16nm (A 2 , 
or A 4 , after aging), 2.33 nm {B 2 ), 2.5 nm (Ci) and 2.9 nm (Fq). The resistivity was measured as a 
function of temperature with various magnetic field ranging from OT to 14 T and for sample A 2 
with various applied currents in the vicinity of the QCP 

The critical temperature was determined by using the extrapolation to zero resistance of the 
tangent at the inflection point of the resistance versus temperature curve. See the values in Table 1 
of the Appendix. 

High magnetic field measurements were performed at LNCMI (Toulouse) on sample Fq, using 
pulsed fields up to 55 T. The ac resistance was measured using four probe technique at a frequency 
of 50 kHz. This was used to perform scaling analysis on sample Fq at the vicinity of the QCP. 

In order to determine the paraconductance contribution, one needs to extract the normal state 
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conductance. For samples A 2 , A 4 and Yq, that were close enough to the insulating transition, so 
that the normal state was recovered under a perpendicular 14 T magnetic field, we identified aN{T) 
with ai 4 T'{T). For the thicker samples B 2 , Ci and Fq, the normal state was not recovered under 
14 T and we extrapolated the normal state from the high-temperature data (typically between 80K 
ad 150K) using the same fitting law that worked for the 14 T resistance of the thinner samples, but 
with different parameters, namely Pn{T) = + aiLx + oq, with Lt = ln(l/T). (For Xq, 

that was measured in situ inside the STM cryostat, a similar extrapolation was used for the normal 
state.) This extrapolation may induce an error on the absolute amplitude of Aa, but not on the 
exponent governing the temperature dependence of the divergent term in (t(T) at T = Tc, since 
no divergence at Tc is present in the normal state conductance nor in its fitting law. 
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APPENDIX 


SAMPLE FABRICATION AND CHARACTERISATION 

The NbN films were grown at Karlsruhe Institute of Miero and Nanoeleetronie Systems, on the 
optieally polished side of sapphire substrates by means of do reaotive magnetron sputtering of a 
pure Nb target in an Ar + N 2 gas mixture as desoribed in Ref. [l251l . The deposition prooess was 
optimized in order to provide the highest transition temperature. 

The thiokness was evaluated during the deposition of NbN at the Karlsruhe Institute of Mi¬ 
cro and Nanoeleetronie Systems, by time of exposure of the substrate and measurement of the 
deposition rate. Figurej^ shows a high-resolution transmission electron microscopy image of a 
2.33 nm-thick NbN film. The atomic rows of the sapphire substrate and of the NbN layers are 
clearly visible. The NbN film consists of nanocrystals of typical size of approximately 2-5 nm, 
linked together by crystalline boundaries. 

Both the chemical analysis performed across the film and the atomic force microscopy (AFM) 
measurements indicated that an oxide layer (presumably Nb 205 ) is present on top of the NbN 
films. This layer appears within a few minutes after exposing the film to air for the first time. The 
thickness of the oxide layer on newly prepared films varies in the range 0.5 to 1 nm. Further oxi¬ 
dation proceeds much slower. The transition temperature of samples with Tp < 10 K was actually 
found to decrease with time and thermal cycling. Simultaneously, the normal-state resistance was 
found to increase. This fact was used experimentally to tune the superconductor/insulator tran¬ 
sition in sample A 4 which corresponds to aged A 2 . We were able to establish that the top oxide 
layer is insulating and not conducting. The dependency of the STM tunneling current versus tip- 
sample distance shows an exponential decay whose exponent is directly linked to the effective 
work function of the sample. For a metallic sample, it is typically around 4eV. When an ultra- 
thin (ie few atomic monolayers thick) insulating barrier on top of a metal is present, the effective 
work function is lowered (see [|48l ). typically by a factor of two at least. We have observed that 
the dependency of the tunneling current shows such a reduced work function, thus supporting the 
existence of an insulating oxide layer. Therefore the existence of a top metallic layer such as NbO 
for instance seems unlikely in our case. The effect of this oxide layer on the superconducting 
properties of ultra thin NbN films is discussed in more details in [|^ . 
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Figure 5. High-resolution transmission electron microscopy (HRTEM) image of a 2.33 nm-thick NbN film 
grown on sapphire substrate. The NbN film consists of nanocrystals of typical size of approximately 2-5 nm, 
linked together by crystalline boundaries. The presence of a ~ 1 nm thick insulating oxide layer (possibly 
Nb20^) on top of the NbN film was inferred from a chemical analysis with scanning TEM measurements 
across the layers and from AEM measurements and is not visible on this image. 

SCANNING TUNNELING SPECTROSCOPY (STS) DATA 

We present in Figj^) and b) the representative dI/dV{V) tunneling spectra measured at 
300 mK and 4.2 K on the area of interest shown in Fig. la of the main text. The blue spectra are 
measured on the lower gap areas (blue areas) of the gap map presented in Fig. lb of the main text, 
while red spectra are measured on larger gap areas of Fig. lb. It is seen that a spatially varying 
pseudogap is measured at 4.2 K, having the characteristic energy scale of the superconducting 
gap measured at low temperature. Figs.[^) and d) present dl/dV conductance maps integrated 
in the energy range of the low-temperature superconducting gaps, at 300 mK (c) and at 4.2 K (d). 
It is clear from this figure that both superconducting correlations and superconducting inhomo¬ 
geneities persist above Tc- In addition, our analysis shows that these two conductance maps c) 
and d) measured over the same topographic area are strongly correlated. This result can be seen 
in Fig. [^) where the cross-correlation map between |^) and |^) is plotted. This proves that the 
inhomogeneities in the energy gap value below Tc and in the pseudogap features above Tq are 


18 



Figure 6. dl/dV tunneling conductance spectra as a function of bias voltage measured by STS on sample 
Xq at a) 300 mK, b) 4.2 K. A spatially varying pseudogap is measured at 4.2 K, having the characteristic 
energy scale of the superconducting gap. The blue spectra are measured on the lower gap areas (blue areas) 
of the gap map shown in Fig. lb of the main text, while red spectra are measured on larger gap areas of 
Fig. lb. c) dl/dV conductance map integrated in the energy range ±0.9 meV (inside the low-temperature 
energy gap region) at 300 mK. The observed inhomogeneities correspond to those seen in the gap map of 
the main text shown in panel Fig. lb. d) Same quantity as in c) but at 4.2 K. This map evidences inho- 
mogenities above Tc in the pseudogap features, e) Cross-correlation map between c) and d) showing that 
the inhomogeneities at 300 mK and at 4.2 K are strongly correlated, f) Radial profile extracted from the 
cross-correlation map e). This analysis allows us to conclude that the length scale Lj characteristic of the 
superconducting inhomogeneities is the same at 300 mK and at 4.2 K. 

strongly spatially correlated. Furthermore, the radial profile shown in Fig.[^, extraeted from a eir- 
cular average of the cross-eorrelation map e), allows us to infer that the eharaeteristie correlation 
length Li of these inhomogeneities at 4.2 K is eomparable to the one seen at 300 mK, of size about 
100 nm. We have also performed STS measurements on sample Xq at higher temperature so well 
in the OD regime (around 7K) whieh show similar results as the one presented in Fig.[^) and Fig. 
1^) at 4.2K. The only differenee between the 7K and 4.2K data is that the speetroseopie features 
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related to the pseudogap are smaller in amplitude at 7K. But the regions where the pseudogap is 
still present (eorresponding to the supergrains) have the same eharaeteristie size and spaeing as at 
300mK or at 4.2K. 


ANALYSIS OF THE PARACONDUCTIVITY : FITTING PARAMETERS FOR THE NORMAL 
STATE RESISTANCE 

For samples A 2 , A 4 and Yq, that were elose enough to the insulating transition, so that the 
normal state was reeovered under a perpendieular 14 T magnetie field, we identified a^iT) with 
o'uTiT). The resistanee of sample Xq was measured in situ inside the STM eryostat up to 30 K, 
therefore the measure under 14T magnetie field eould not be performed, and we extrapolated the 
normal state. For the thieker samples, B 2 , Ci and Fq, supereonduetivity eould not be eompletely 
suppressed by a magnetie field of 14 T. In sueh eases, the normal-state resistanee was fitted by an 
expression that worked well for the thinner samples under 14T, 


Pn{T) — + UiLt + a^, (4) 

with Lt = ln(l/T). We determined the parameters 09 , 1,3 fitting the high-temperature resis¬ 
tanee, typieally in the range 80— 150 K (exeept for the sample Xq, for whieh the range was foreedly 
smaller, 10—30K). The fitting parameters are found in Tab.|^ 

ANALYSIS OF THE PARACONDUCTIVITY : POSSIBLE OBSERVATION OF A 0D-2D CROS¬ 
SOVER 

In Fig.|^we plot the pararesistanee Aa~^ as a funetion of temperature (and not redueed tem¬ 
perature e) for sample Xq. Sueh plot (similar the right inset of Fig. 3b in the Artiele for sample 
A 2 ) allows to visualize both the OD AL (quadratie) regime and the 2D (linear) regime eloser to 
Tc- Sinee the two regimes are eharaeterized by different “fluetuational” eritieal temperatures, they 
eannot be both visualized with a single ehoiee of Tc in the redueed temperature e. Sueh 0D-2D 
erossover was indeed observed in all thinner samples A 2 , A 4 , Xq and Yq. 
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sample name 

thickness 

(nm) 

Tc 

(K) 

i?sg(300K) 

(kD/sq) 

(K) 

(K) 

Fo 

2.90 

9.0 

1390 

8.6 

- 

Cl 

2.50 

9.4 

2370 

9.5 

- 

B 2 

2.33 

7.1 

2450 

7.1 

- 

^2 

2.16 

4.5 

3150 

4.9 

4.4 

A 4 (aged A 2 ) 

2.16 

2.4 

4570 

3.0 

2.9 

Xo 

2.14 

3.8 

3800 

3.9 

3.6 

Yo 

2.10 

4.3 

2350 

4.9 

4.3 


Table I. Sample nominal thickness, critical temperature Tc (defined as the extrapolation to zero resistance 
of the tangent at the inflection point of the resistance curve), resistance per square Rgq at room temperature 
for the different samples studied, and “fluctuational” critical temperature for the 2D (T^^) and OD (T^^) 
AL regime. 


sample name 

ao (kD/sq) 

ai (kD/sq) 

as (kD/sq) 


9299 

1038 

1.370 

B 2 

5804 

275.3 

10.53 

Cl 

3529 

32.32 

6.094 

Fo 

2212 

4.512 

4.381 


Table II. Fitting parameters of the normal state resistance according to the expression pn{T) = + 

“2-^t + + 0 , 0 , with Lt = ln(l/r). Since the inclusion of the L|, term did not improve significantly 

the fits, to reduce the number of fitting parameters, we fixed 02 = 0. 

STUDY OF THE QUANTUM CRITICAL REGION 

In Fig.are displayed the isothermal resistanee versus magnetie field curves for different tem¬ 
peratures. The curves for temperatures between 1.8 K and 2.6 K exhibit a common crossing point, 
at He = 10.7 T and Rc = 11 kf2, corresponding to the quantum critical point. For these measure¬ 
ments, the current and voltage values across the sample were very low and the curves were found 
to be insensitive to amplitude variations of the current bias. In Fig.are displayed the isopoten¬ 
tial resistance versus magnetic field curves for different values of the voltage applied across the 
sample, corresponding to different values of the electric field. The measurements were performed 
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Temperature (K) 


Figure 7. Pararesistance Aa~^ for sample Xq as a function of temperature, illustrating the OD AL be¬ 
havior, where Acr“^ oc (T — and the crossover to the 2D AL regime close enough to Tc, where 

Af7“^ oc T — . We point out that the “fluctuational” critical temperatures and need not be the 

same, as they enter the expressions of fluctuation contributions living on different length scales. 



Magnetic Field (T) 

Figure 8. Resistance versus magnetic field data for temperature between 1.8 K (lower curve at low fields) 
and 2.6 K (upper curve at low fields). 


at 1.9 K. As for the isothermal eurves, the isopotential eurves are found to eross in a single point. 
These eurves were then used for sealing analysis as deseribed in the Artiele. 
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Figure 9. Resistance as a function of magnetic field at different electrical fields across the sample at a 
temperature of 1.9 K. 

THEORETICAL ANALYSIS OF PARACONDUCTIVITY 

Zero dimensional Aslamazov-Larkin (AL) fluctuations 

To discuss AL paraconductivity, we assume overall isotropy (at least in the average) and start 
from the expression [|2^ |4^ 



(5) 


that reproduces the standard AL result in D dimensions with 7 q = being the 

relaxation rate for each mode q. is the coherence length and e = log(T/Tc), Tc being the 
critical temperature seen by the fluctuating Cooper pairs (which may differ from the actual critical 
temperature, in the presence of crossover phenomena). The above formula is rather general fISOl . 
although the expression of e may change depending on the microscopic theory, and the expression 
e = a{T — Tc)/Tc, with a dimensionless prefactor a not necessarily equal to one, might be used. 
However, for definiteness, we adopt henceforth the standard definition of e. 

Since 7 q only depends on q, the angular integration factors out and one obtains 



( 6 ) 


where = 271^^'^/T{D/2) is the surface of the unitary sphere in D dimensions and r( 2 ;) is the 
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Euler gamma function. Exploiting the expression of 7 ^ one obtains 


Sao = 


ire'^QD^o 


dq 


q 


D+l 


hD{27i)^ J, (e + CoV) 


\3 • 


(7) 


When D ^ 0, Qd inherits the singularity of r(Zi)/2) ^ 2/D, so that Qd 


D, and 




( 8 ) 


4^62 ’ 

i.e., there exists a well defined limit of the AE paraconductivity as D —0, that scales with e“^, and 
qualitatively agrees with the behavior observed in the thinnest NbN films. This result is converted 
into the measured conductance per square, dividing it by the square of a suitable length scale I, 
related with the size of the fluctuating (nano)domains, 


(5 a. 


measured 


/- \ 2 2 

^o\ Tre 


(9) 


I I Ahe^' 

However, although the supergrains may behave as quantum dots, closer to Tc the system is two- 
dimensional, the supergrains are coupled, and the paraconductive fluctuations must cross over to 
D = 2. The formal limit does not allow for a description of this crossover and the identification 
of the relevant length scales. A crucial remark of Ref. P9l is that the AE theory can be cast in an 
even more general form, introducing a suitable weighted density of states, 

d^q 




j (27r)^ 

which allows to write the paraconductivity in the form 


(Vq7q) 5(7 -7q), 


( 10 ) 


5(Td = 


7re 

4^ 


d7 


(e + 7)3 


( 11 ) 


The standard AE result is recovered with 


yielding 


■A/d(7) = 


dan = 


dq 


(27r)^ Jo 


D 


d7 


7 


D/2 


( 12 ) 


(13) 


2hD{2n)D -'(e + ^)3- 

We argue that the observed 0D-2D crossover stems from a different form of A/"d(7), that results 
from the supergrains (endowed with an internal structure) being connected in a 2D network. In the 
following we thus set D = 2. The standard 2D AE result is found for M 2 ( 7 ) = 7 /tt, corresponding 
to diffusion of fluctuating Cooper pairs in two dimensions. If instead A 2 ( 7 ) = 7 / 7 r, for 7 < 7 , 
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and A 2 ( 7 ) = 7 / 7 r, for 7 > 7 , which describes a slowing-down of the diffusion of fluctuating 
Cooper pairs above a threshold 7 (e.g., corresponding to trapping of the fluctuations within the 
supergrains), the paraconductivity is 


5a 


7 

\Qh e(e + 7 ) ’ 


(14) 


which recovers the 2D behavior for e 7 , while for e 3> 7 

16 ^ 62 ’ 


(15) 


mimicking the OD behavior at higher temperature. A comparison with the formal extrapolation 
of AL fluctuations for D —)■ 0 yields 7 = 47r(.^o/0^- ^ matter of fact, the 0D-2D crossover is 

clearly visible for sample Xq (d = 2.14 nm) for which a much greater number of points was taken 
in the vicinity of the transition, as well as for the sample A 2 (d = 2.16 nm), but also for A 4 and Yq. 
Eventually, at temperatures below the 2D AL regime, the system may be governed by percolation 
physics or, alternatively, by Berezinski-Kosterlitz-Thouless behavior. 

A rough estimate gives e ~ where ^ is the temperature dependent correlation length 

that diverges at the transition. Then, the 0D-2D crossover takes place at a temperature T such that 
^{T) Ri 'CoV7- A comparison with the formal OD limit of the AL theory allows us to identify 
'j ^ Att hence the size of the supergrains I ^ 
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